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We derive in the present work topological photonic states purely based on silicon, a conventional 
dielectric material, by deforming a honeycomb lattice of silicon cylinders into a triangular lattice of 
cylinder hexagons. The photonic topology is associated with a pseudo time reversal (TR) symmetry 
constituted by the TR symmetry respected in general by the Maxwell equations and the Cq crystal 
symmetry upon design, which renders the Kramers doubling in the present photonic system with the 
role of pseudo spin played by the circular polarization of magnetic field in the transverse magnetic 
mode. We solve Maxwell equations, and demonstrate new photonic topology by revealing pseudo 
spin-resolved Berry curvatures of photonic bands and helical edge states characterized by Poynting 
vectors. 
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Introduction .—The discovery of quantum Hall effect 
(QHE) opened a new chapter of condensed matter 
physics with topology as the central concept [inni. 
Topological states are not only interesting in an academic 
point of view, but also expected to yield significant im¬ 
pacts to applications because robust surface (or edge) 
states protected by bulk topology provide new possibil¬ 
ities for spintronics and quantum computation [iM3- 
However, topological matters confirmed so far are still 
limited in number, and most of them exhibit topological 
properties only at very low temperatures, which hinders 
their better understanding and manipulation indispens¬ 
able for practical applications. 

Photonic crystals are optical analogues of solids with 
lattice of atoms replaced by medium of periodic electric 
permittivity and/or magnetic permeability 0 . Meta¬ 
materials are designed to generate electromagnetic (EM) 
properties such as negative index, magnetic lens, and so 
on, which are not available in nature [l9|. Recently it 
has been recognized that topological states characterized 
by unique edge propagations of EM wave can be real¬ 
ized in photonic crystals based on gyromagnetic mate¬ 
rials under external magnetic field, bi-anisoctroic meta¬ 
materials with coupled electric and magnetic fields where 
bi-anisotropy acts as effective spin-orbit coupling, and 
coupled resonator optical waveguides (CROWs) |2QI431|. 


(for a review see j32|). 

In the present work, we propose a two-dimensional 
(2D) photonic crystal purely made of silicon, a conven¬ 
tional dielectric material. We notice that a honeycomb 
lattice is equivalent to a triangular lattice of hexagonal 
clusters composed by six neighboring sites, and that, tak¬ 
ing this larger hexagonal unit cell instead of the primitive 
rhombic unit cell of two sites (see Eig. Hl^a)), the Dirac 
cones at K and K’ points in the first Brillouin zone of 
honeycomb lattice are folded to doubly degenerate Dirac 
cones at T point. It is then intriguing to observe that 
at the r point there are two 2D irreducible representa¬ 



FIG. 1. (a) Schematics of a triangular photonic crystal of 

’’artihcial atoms” composed by six cylinders of dielectric ma¬ 
terial with dielectric constant ea along the z direction embed¬ 
ded in the environment of dielectric constant ca • Red dashed 
rhombus and hexagon are primitive cells of honeycomb and 
triangular lattices respectively. Solid black hexagon labels an 
artificial atom, while dashed black one marks interstitial re¬ 
gions among atoms. di and a 2 are lattice vectors, and ao 
is the lattice constant of triangular lattice. Right panel: en¬ 
larged view of hexagonal cluster with R the length of hexagon 
edge and d the diameter of cylinders, (b) Pseudo spin states 
of the present photonic system associated with the left- and 
right-hand circular polarizations of magnetic fields Hx ± iHy 
of the TM mode. 


tions in Cq symmetry group associated with odd and even 
parity respective to spatial inversion operation. Based on 
these properties, we propose to open a topologically non¬ 
trivial band gap by deforming the honeycomb lattice in 
the way keeping hexagonal clusters and preserving the 
Cq rotation symmetry (see Eig. [T]^a)). Explicitly, we re¬ 
veal by solving the Maxwell equations that harmonic EM 
modes hosted by the hexagonal cluster, working as ”ar- 
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tificial atom” in the present scheme, exhibit electronic 
orbital-like s-, p-, and d-wave shapes and form photonic 
bands. We clarify that there is a pseudo time reversal 
(TR) symmetry constituted by the TR symmetry re¬ 
spected by the Maxwell equations and the Cq crystal 
symmetry upon design, which behaves in the same way 
as TR symmetry in electronic systems and renders the 
Kramers doubling in the present photonic system. This 
intimately gives the one-to-one correspondence between 
the left- and right-hand circular polarizations of magnetic 
field in the transverse magnetic (TM) mode and the up- 
and down-spin states of electrons. Evaluating the Berry 
curvatures of photonic bands and the edge states for fi¬ 
nite systems, we demonstrate the emergence of topolog¬ 
ical phases. With the simple design free of requirement 
on any external field and gyromagnetic or bi-anisotropic 
materials, the present topological photonic states purely 
based on silicon are expected very promising for future 
applications. 

Artificial atom and pseudo spin. —Let us consider har¬ 
monic TM modes of EM wave, namely those of finite 
in-plane Hy^ out-of-plane Ez components and with 
others being zero, in a dielectric medium (for coordinates 
see Eig.[Ha)). Eor simplicity, the real electric permittiv¬ 
ity is taken frequency independent in the regime under 
consideration. The master equation for a harmonic mode 
of frequency uj is then derived from the Maxwell equa¬ 
tions [33! 


V X Vx 

£(r) 





( 1 ) 


with 5(r) the position-dependent permittivity and c the 
speed of light. The transverse components of magnetic 
field are given by the Earaday relation H = —i/(/io<^) V x 
E, where the magnetic permeability /tq is presumed as 
that of vacuum. The Bloch theorem applies for the 
present system when ^(r) is periodic in space as shown 
schematically in Eig. [IJa). It is reminded however that 
the master equation describes the EM waves instead 
of electrons carrying on the spin degree of freedom, with 
the most prominent difference lying at the response upon 
TR operation. 

Our photonic crystal is made of cylinders of dielectric 
material with relative permittivity Sd (such as silicon) 
parallel to the 2; axis embedded in a dielectric background 
of £a (such as air) as schematically shown in Eig. [U^a). 
Six dielectric cylinders of diameter d form a hexagon with 
edge length R. Our discussions below are for states uni¬ 
form along the ^ axis, which reduces the problem to 2D. 
We start from a honeycomb lattice of dielectric cylinders, 
and deform it in the way keeping hexagonal clusters com¬ 
posed by six neighboring cylinders and the Cq symmetry. 
Now the alignment of dielectric cylinders is more conve¬ 
nient to be considered a triangular lattice of hexagonal 
’’artificial atoms”. There are two 2D irreducible repre¬ 
sentations in the Cq symmetry group associated with the 



FIG. 2. (a) Electric field Ez of lYiepx/py and dxy jd^^_y 2 pho¬ 
tonic orbitals hosted by the artificial atom, (b) Two pseudo 
spin states corresponding to the left- and right-hand circular 
polarizations of magnetic fields Hx ± iHy associated with Ez 
fields of positive and negative angular momenta: p± and d±. 


triangular lattice: E' and E" with basis functions x/y 
and xyj{x^ — ^^), corresponding to odd and even spa¬ 
tial parities respectively [3^. As can be seen from Ez 
fields in Eig. [2](a) obtained by solving Eq. m, the artifi¬ 
cial atoms cany px/py and dxyldx 2 -y 2 orbitals, with the 
same symmetry as those of electronic orbitals of conven¬ 
tional atoms in solids. 

We now examine matrix representations of 7r/3 rota¬ 
tion and its combinations for the basis functions of PxjPy 
and dxyldx 2 -y 2 . Since PxjPy behave in the same way as 
x/y^ it is easy to see 


It is noticed that U = [De'{Cq) + De'{Cq)]/VS = —iCFy 
with De'{Cq) = D‘^^,(Cq) is associated with the 7r/2 ro¬ 
tation of PxtPy {cFy being the Pauli matrix). Therefore, 
U‘^{Px^Py)^ = —{Px^Py)^^ which is consistent with the 
odd parity of PxjPy with respect to spatial inversion. 
Similarly, one has 

= )•<='> 

which, in the matrix form, is same as De>{Cq) be¬ 
cause the basis functions are now bilinear of x/y. 
It is then straightforward to check that [De"{Cq) — 
De"{CI )\/= u is associated with a 7r/4 rota¬ 
tion of dxy/dx^-y-^^ which yields U‘^{dx‘ 2 -y 2 ^dxy)^ = 

(dj,2_y2 , ) 

We compose the antiunitary operator T = UK, where 
K is the complex conjugate operator associated with the 
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FIG. 3. Dispersion relations of TM mode for the 2D pho¬ 
tonic crystals with = 11.7, sa = 1 and d = 2R/3 for (a) 
ao/R = 3.16 (Inset: Brillouin zone of triangular lattice), (b) 
ao/R = 3 and (c) ao/R = 2.78. Blue and red are for d± 
and p± bands respectively, and rainbow for hybridization be¬ 
tween them. The case of ao/R = 3 corresponds exactly to 
the honeycomb lattice of individual cylinders. Diameter of 
individual cylinders is kept the same for all calculations since 
its variation influences photonic bands slightly in quantitative 
manners. 


TR operation respected by Maxwell systems in general. 
Since = —1 is guaranteed by 77^ = —1, T can be 
taken as a pseudo TR operator which provides Kramers 
doubling in the same way of electronic systems. It is clear 
that the crystal symmetry plays an important role in this 
pseudo TR symmetry [35|. 

The two pseudo spin states are given by 

P± = {Px ± iPy)/ \/2; d± = {dx2_y2 ± id^y)/ \/2, (4) 

which are related to the above basis functions by uni¬ 
tary transformation. Namely the pseudo spin up and 
down correspond to the wave functions of fields of 
positive and negative angular momenta, or equivalently 
the left- and right-hand circular polarizations of in-plane 
magnetic fields defined by H^EiUy as shown in Figs.^b) 
and[2fb). 

Pseudo spins discussed so far in photonic systems in¬ 
clude bonding/antibonding states of electric and mag¬ 
netic fields [2J, [23, left-hand/right-hand circular polar¬ 
izations of EM waves [28|, and clockwise/anticlockwise 
circulations of light in CROWs (2^.1^. The first two im¬ 
plementations involve sophisticated metamaterial struc¬ 
tures and thus require delicate fabrications. While the 
last one only uses silicon fibers, and thus easy to prepare, 
a break in one of coupling loops at perimeter induces a 
’’magnetic impurity” which destroys the topological he¬ 
lical edge states. 

Photonic bands .—Now we calculate the photonic band 
dispersions described by the master equation m impos¬ 
ing periodic boundary conditions along unit vectors ai 
and 02 given in Fig. [IJa). As shown in Fig. [H double 
degeneracy in the band dispersions appears at F point, 
which can be identified as p± and d± states, consistent 
with the symmetry consideration. For large lattice con¬ 
stant ao, the photonic band below (above) the gap is oc- 



FIG. 4. Real-space distributions of the time-averaged Poynt- 
ing vector associated with the pseudo spin-down state at F 
point below the photonic gap: (a) ao/R = 3.16 in the trivial 
regime and (b) ao/R = 2.78 in the topological regime. 


cupied by p± {d±) states (see Fig.[3l^a) for ao/R = 3.16). 

Reducing the lattice constant to ao/R = 3, the p and d 
states become degenerate at F point, and two Dirac cones 
appear as shown in Fig.^b). This is because that at this 
lattice constant the system is equivalent to honeycomb 
lattice of individual cylinders, and the doubly degenerate 
Dirac cones are nothing but those at K and K’ point 
in the Brillouin zone of honeycomb lattice based on the 
primitive rhombic unit cell of two sites 0. 

When the lattice constant is further reduced, a global 
photonic band gap is reopened as shown in Fig. [3fc) for 
ao/R = 2.78. Now the Ez field at low-(high-)frequency 
side of the band gap exhibits d±{p±) characters around 
the F point, opposite to the order away from the F point. 
Namely a band inversion takes place upon reducing the 
lattice constant in the present photonic lattice. Quanti¬ 
tatively, the band gap is Auj = 13.8 THz at uj = 138.75 
THz with ao = 1 /im, with all the quantities scaling with 
the lattice constant. 

In order to see what happens in the system around 
the band inversion, we check the real-space distribution 
of Poynting vector S = Re[R x R*]/2 averaged over a 
period r = 27r/ uj, which describes the energy flow in 
the present EM system. It is found that the Poynting 
vector is circling around individual atoms as shown in 
Eig. m^a) for ao/R = 3.16, with the chirality of Poynting 
vector corresponding to the pseudo spin (Poynting vector 
with pseudo spin-up is not shown explicitly). The EM 
energy flows around individual atoms, characterizing a 
conventional ’’insulating” state. At ao/R = 2.78, namely 
after the band inversion, the Poynting vectors are much 
enhanced in interstitial regimes as shown in Eig.lU^b). It 
is in a sharp contrast to the case in Eig.lU^a), and implies 
an unconventional insulating state. 

Topological edge state .—We also consider a ribbon of 
photonic crystal after band inversion by cladding its two 
edges in terms of two photonic crystals with trivial band 
gap (namely before band inversion) at the same fre¬ 
quency window, which prevents possible edge states from 
leaking into free space. It should be kept in mind that, 
since the cluster of six cylinders is the basic block of the 
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FIG. 5. (a) Dispersion relation of a ribbon-shaped topological 
photonic crystal, which is infinite in one direction and of 45 
and 6 artificial atoms for the topological and trivial regions 
respectively in the other direction. Right panel: enlarged 
view of (a). Red curves are for topological edge states. A and 
B are two points where Ez field are shown in (b). Parame¬ 
ters used here are same as those in Fig. [3] except ao/R = 2.9 
and ao/R = 3.125 in topological and trivial regions respec¬ 
tively. (b) Real-space distributions of Ez fields at points A 
and B. Right panels: time-averaged Poynting vectors S over 
a period, (c) Photonic crystal of height h with two horizon¬ 
tal gold plates placed at two ends symmetrically. Hexagons 
in fade white color are primitive cells for the triangular lat¬ 
tice. (d) Distribution of energy-density carried by Ez held 
uezM — /2 in a topological photonic crystal 

cladded by a trivial one with h = 10 mm, d = 2.4 mm and 
ao = 10 mm, and R = 3.65 mm and 3 mm in topological and 
trivial regions respectively. Frequency of harmonic line source 
is cj = 13.47 GHz within the photonic band gap. 


present design, it should not be destroyed for any mean¬ 
ingful discussion. As displayed in Fig.[5fa), there appears 
edge states as indicated by the double degenerated red 
curves, with a small gap (not noticeable in this scale) 
opened at F point due to the suppression of Cq sym¬ 
metry and then pseudo TR symmetry at the interface. 
Checking real-space distribution of Ez field at typical mo¬ 
menta around F point (A and B in the enlarged vision 
of Fig. [5ja) with kx = ±0.04|^), we find that the in-gap 
states locate at the edge and decay exponentially into 
bulk as displayed in Fig. [5fb) (two other states are local¬ 


ized at the other ribbon edge and not shown explicitly). 
As shown in the right insets of Fig. [Ifb), the Poynting 
vectors exhibit a nonzero downward/upward EM energy 
flow for the pseudo spin-down/-up state even averaged 
over time. This indicates unambiguously counter propa¬ 
gations of EM energy at the sample edge associated with 
the two pseudo spin states, the hallmark of a quantum 
spin Hall effect (QSHE) state [ 1 , Distributions of 
Poynting vectors of the bulk bands in Eig. for the 
ribbon system are similar to those in Eig. lU^b) for the 
infinite system. QHE has been described by the cyclic 
motions of electrons under strong external magnetic field 
in a quasi classic picture of electronic wave functions . 
It is noticed that the Poynting vector is a physical quan¬ 
tity in EM systems, and therefore the distributions shown 
in Eigs. |l]and[5fb) can be observed in experiments. The 
photonic QSHE in the present system can also be con¬ 
firmed by evaluating the Z2 topology index based on a 
k • p model around F point. Although Dirac dispersions 
in photonic systems were discussed previously in both 
square and triangular lattices j37l-[39| , possible nontrivial 
topology was not addressed. 

For experimental implementation of the present topo¬ 
logical state, the finite height of silicon rods along z di¬ 
rection has to be taken into account. We consider ex¬ 
plicitly a square sample of topological photonic crystal 
sandwiched by two horizontal gold plates separated by 
h (see Fig. [5fc)), with the height h chosen to prevent 
photonic bands with nonzero kz from falling into the 
topological band gap. The size of topological sample is 
40ai X 20(ai + a2) with all four edges cladded by a trivial 
photonic crystal. A harmonic line source E = Eqc^^'^z 
is placed parallel to silicon rods to inject EM wave at 
the interface with the frequency in the topological band 
gap. We simulate the system by solving time-dependent 
Maxwell equations j4l|, |42|. Since any harmonic source 
preserves TR symmetry respected by Maxwell equations, 
the system exhibits helical topological edge states as 
shown in Eig. [5fd). 

In conclusion, we derive a two-dimensional photonic 
crystal with nontrivial topology purely based on silicon, 
a conventional dielectric material, simply by deforming 
honeycomb lattice of silicon cylinders. A pseudo time re¬ 
versal symmetry is constructed in terms of the time rever¬ 
sal symmetry respected by the Maxwell equation in gen¬ 
eral and the Cq crystal symmetry upon design, which en¬ 
ables the Kramers doubling with the role of pseudo spin 
played by the polarization of magnetic field of transverse 
magnetic mode. The present topological photonic crystal 
with simple design backed up by the symmetry consider¬ 
ation can be fabricated relatively easy as compared with 
other proposals, and is expected to leave impacts to the 
topological physics and related materials science. 
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